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Abstract 

The Riesz distribution for real normed division algebras is derived in this work. 
Then two versions of these distributions are proposed and some of their properties 
are studied. 



1 Introduction 



In the context of the Statistics and related fields, the analysis on symmetric cones has been 
developed in depth through the last seven decades. As usual, the studies were focused on 
the real case, later the complex case appeared and in the recent years, the quaternionic and 
octonionic numbers were included in the analysis. This analysis has been notorious in the 
development of certain statistical areas as random matrix distributions, but also in some 
mathematical topics including special polynomials matrix arguments as zonal, Jack, Davis' 
invariants, Laguerre, Hermite, Hayakawa; and the associated hypergeometric type series. 
These mathematical tools have been implemented for example in the characterisation of 
non-central distributi ons, s t atistical theory of s h ape, et c . Among many au t hors, we can 
highli g ht the work s of Herz ( 1955 ). James ( 196D . James! (1964), Constantind ( 1963 ). Davis! 
(|l98dh . iMuirheadl (Il982h . iDiaz-Garcial faOlK and the references therein. Of course, the 
interest in the analysis of symmetric cones does not belong exclusively to statistics; the 
harmonic analysis, for example, is an special case of the real symmetric cone and this 
cone play a fundamental role in number theory. Another particular case that has been 
studied in detail, specially from the point of view of the wave equation, is the Lorentz cone, 
Faraut and Koranvil (| 1994T) and the references therein. 

Although during the 60's and 70's were obtained important results in the general anal- 
ysis of the symmetric cones, the past 20 years have reached a substantial development. 
Essentially, these advances have been archived through two approaches based on the the- 
ory of Jordan algebra and the real normed division algebras. A basi c source of the gen- 
eral th eory of symmetric cones under Jordan algebras can be found in iFaraut and Koranvi 
( 1994J ); and specifically, some works in the co ntext of distrib u tion theory in symmetri c 
cones based on Jordan al gebra s are provided in Massaml ( 19941 ). Casalis and Letac ( 1996 ). 
Hassairi and Laiinl ( 200ll ). and lHassairi et al.l(|2005l ). and the references therein. In the field 
of spherical functions (Ja ck polynomials including James' zonal polynomials) we can mention 
the work of Sawyer ( 1997t ) . Parallel results on distribution theory based on real no rmed divi 



sion algebras have been also developed in statistics and random matrix theory, see lForrester 
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(2009|), Di'az-Garci'a and Gutierrez- Jaimez ( 2011 ). Diaz-Garcia and Gutierrez- Jaimej ( 2012a ). 
among others. A detailed study of functions with matrix arguments on symmetric cones, 
such as Jack a nd Davis' invariants polynomials and hypergeometric types series , among oth- 
ers, appear in Gross and Richarda ( 19871) , iDiaz-Garcia and Gutierrez- Jaimezl (|2012bh and 



Diaz-Garcia and Gutierrez- Jaimea (|2009j), for real normed division algebras. 



Under the approach based on the theory of Jordan algebras, a famil y of distributions on 
symm etric cones, termed the Riesz distributions, was first introduced bv lHassairi and Laimi 
(|200lf ) under the name of Riesz natural exponen tial family (Riesz NEF); it was based on 
a specia l case of the so-called Riesz measure from Faraut and Korarivil ( 1994 . p. 137), going 
back to Riesz (1949). This Riesz distribution generalises the matrix multivariate gamma 
and Wishart distributions, containing them as particular cases. 

This article introduces the Riesz distribution for real normed division algebras. Section 
[2] revises some definitions and notations on real normed division algebras, also, two defini- 
tions of the generalised gamma funcion on symmetric cones are given; in addition, several 
Jacobians with respect to Lebesgue measure for real normed division algebras are proposed. 
Finally, two Riesz distributions are obtained and under both definitions, the corresponding 
characteristic function and the eigenvalue distribution are derived in section [3] 



2 Preliminary results 



A detailed discussion of r eal normed division algebras may be found in iBaed (|2002l) and 
Gross and Richards! (|l987h . For convenience, we shall introduce some notation, although in 
general we adhere to standard notation forms. 

For our purposes, a vector space is always a finite-dimensional module over the field 
of real numbers. An algebra $ is a vector space that is equipped with a bilinear map 
m : 5x5 -> 3 termed multiplication and a nonzero element 16 5 termed the unit such 
that m(l, a) — m(a, 1) = a. As usual, we abbreviate m(a, b) — ab as ab. We do not assume 
5 associative. Given an algebra, we freely think of real numbers as elements of this algebra 
via the map u) i— > uil. 

An algebra 5 is a division algebra if given a, b € 5 with ab = 0, then either a = or 
b = 0. Equivalently, 5 is a division algebra if the operation of left and right multiplications 
by any nonzero element is invertible. A normed division algebra is an algebra 5 that is 
also a normed vector space with ||a6|| = ||a||||6||. This implies that 5 is a division algebra 
and that 1 1 1 1 = 1 . 

There are exactly four normed division a lgebra s : real numbers (3?), complex numbers 
(£), quaternions and octonions (D), see iBaezI (|2002l) . We take into account that 5ft, 
£, S) and O are the only normed division algebras; moreover, they are the only alternative 
division algebras, and all div ision algebras have a real dimension of 1,2,4 or 8, which is 
denoted by (3, see IBaezI ( 20021 Theorems 1, 2 and 3). In othe r br anches of ma thematics, the 
parameters a = 2/(3 and t = (3/4 are used, see lSawverl (|1997h and lKabd (|1984l ). respectively. 

Let „ be the set of all m x n matrices of rank m < n over 5 with m distinct positive 
singular values, where 5 denotes a real finite- dimensional normed division algebra. Let 
groxn ^ e se t of all m x n matrices over 5- The dimension of $ mxn over 5ft is (3mn. Let 
A G g m x« 5 tb en A* = A T denotes the usual conjugate transpose. 

Table Q] sets out the equivalence between the same concepts in the four normed division 
algebras. 

We denote by 6^ the real vector space of all S s 3™ xm such that S = S*. In addition, 
let be the cone of positive definite matrices S G g;™x m _ Thus, consist of all matrices 
S = X*X, with X e £^ n ; then ^ is an open subset of 6^. 

Let be the diagonal subgroup of £^ m consisting of all D e ^ mxm ) D = diag(di, . . . , d„ 
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Table 1: Notatii 



Real 



Complex 



Quaternion 



Octonion 



Generic 
notation 



Semi-orthogonal Semi-unitary Semi-symplcctic 



Orthogonal 
Symmetric 



Unitary 
Hcrmitian 



Symplcctic 

Quaternion 
hcrmitian 



Scmi-cxccptional 
type 
Exceptional 
type 
Octonion 
hcrmitian 



el 



For any matrix X 6 ^™ xm j rfX denotes the matrix of differentials (dxij). Finally, 



1 & 13 D 13 



or V^ n , see 



we define the measure or volume eleme nt (dX) when X 6 $ 
Diaz-Garcia and Gutierrez- Jaimea (|201ll ). 

If X g g mxn then (dX) (the Lebesgue measure in J mxn ) denotes the exterior product 
of the f5mn functionally independent variables 



(dX) = A A da;^ where dxij = f\ dx 



(fc) 



fc=i 



If S e ©^ (or S 6 is a upper triangular matrix) then (dS) (the Lebesgue 

3 

m(m — l)/3/2 + m functionally independent variables, 



measure in ©^ or in T^(m)) denotes the exterior product of the exterior product of the 



(dS) = f\ dsn f\ f\ ds 



(fe) 



i<j k=l 



Observe, that for the Lebesgue measure (dS) defined thus, it is required that S £ that 
is, S must be a non singular Hcrmitian matrix (Hcrmitian definite positive matrix). 

If A G then (dA) (the Le gesgue measure in D(^) denotes the exterior product of the 
/3m functionally independent variables 



(dA) = f\ f\d\ 



(fe) 



i=i k=i 

If Hi e V£ |fl then 

(HJdHi) = A A h J dh - 

where H = (HJ|H|)* = (hi, . . . , h m |h m+ i, . . . , h„)* e il^(n). It can be proved that this 
differential form does not depend on the choice of the H2 matrix. When n = 1; V^ n x defines 
the unit sphere in $ m . This is, of course, an (to — l)/3- dimensional surface in $ m . When 
n = to and denoting Hi by H, (HdH*) is termed the Haar measure on il^(m). 

The multivariate Gamma function for the space ©^ denotes as [a] , is defined by 

r&[a] = / etr{-A}|A| -< ff, - 1 > /, / a - 1 (<iA) 

m 

= jr^-iW 4 r[a - (»- 1)0/2], 
»=i 

where etr(-) = exp ( tr(-)) | • | denotes the determinant and Re(a) > (to — l)/?/2, see 
Gross and Richards! (|l987l ). This can be obtained as a particular case of the generalised 
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gamma function of weight re for the space 6^ with re = (fci, k 2 , ■ • • , k m ), k\ > k 2 > • • • > 
fc m > 0, taking re = (0, 0, . . . , 0) and which for Re(a) > (m — l)/3/2 — fe m is defined by, see 



Gross and Richards! |l987t ). 



a, re = 



/ etr^AIIAp-^-^/^VtAXdA) 

m 

^m(m-l)/3/4 JJ p[ fl + fe . _ (• _ 1 ) /3 / 2 ] 



a , 



where for Ae6 



/3 



g„(A) = |A m | fe ™ J] |A, 



fci -fci 



(1) 



(2) 



(3) 



with A „ = (ars), r, s = 1 , 2, p = 1, 2, ...,m is termed the highest weight vec- 

tor, see iGross and Richards! (|19870 . In other branches of mathematics the highest weight 
vector q K (A) is also term ed t he generalised power of A and is denoted as A K (A), see 



Faraut and Koranvil (1994) and Ha ssairi and Laim i (2001). 



Some additional properties of <j K (A), which are immediate consequences of the definition 
of ? K (A) and the following property 1, are: 

1. if A i, ... , A m , are the eigenvalues of A, then 



q K (A) = f[X, 



k, 



(4) 



i=l 



q K (A- 1 )= q - 1 (A)=q_ K (A), 

3. if K = (p, . . . ,p), then 

<Z«(A) = |A|* 
in particular if p = 0, then q K (A) = 1. 

4. if r = (h, t 2 , . ■ ■ , t m ), tx > t 2 > ■ ■ ■ > t m > 0, then 

q K+T (A) = q K (A)q T (A), 
in particular if r = (p,p, ■ ■ • ,p), then 

q K+T (A) = q K+p (A) = \A\Pq K (A). 

5. Finally, for B G S mxm such that C = B*B G ©f„, 

9k(BAB*) = <2«(C) 9k (A) 

and 

<7k(B~ 1 AB* _1 ) = (<z K (C))-V(A). 



(5) 
(6) 

(7) 
(8) 

(9) 
(10) 



Remark 2.1. Let V(&^ n ) denote the algebra of all polynomial functions on &^ nl and 
Vk{&m) the subspace of homogeneous polynomials of degree k and let V K {&^ n ) be an irre- 
ducible subspace of V{&^ n ) such that 

Note t hat q K is a homogeneous polynomial of degree k, moreover q K G 7 ?K (@^ 1 ), see Gross and Richards] 
(|l987h . 
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In @, [a]% denotes the generalised Pochhammer symbol of weight n, defined as 



[af K =f[(a-(i- 1)0/2)^ = ^^, 

i= i r »H 

where Re (a) > (rn—X)(3/2—k m and (a), = a(a+l) • • ■ (a+i— 1), is the standard Pochhammer 
symbol. 

A variant of the generalised gamma function of weight k is obtained from Khatri ( 1966h 
and is defined as 



r£[a,-«] = / etxi-AyW-^-Wfi-^A-^dA) 

m 

= TrTO^-iW^JJrfa- fci — (* — 1)j9/2] 

(-i)*r&[a] 



(ii) 



[-a + (ro - 1)0/2 + 1]£ 



(12) 



where Re(a) > (m - l)/3/2 + fc x . 

Now, we show three Jacobians in terms of t he /3 parameter, which are proposed as ex - 
tensions of real, complex or quaternion cases, see lDi'az-Garci'a and Gutierrez- Jaimea (|2011l ). 

Lemma 2.1. Let X and Y G matrices of functionally independent variables, and let 
Y = AXA* + C, where A £ CJL and C € are matrices of constants. Then 

' III. I It 'III v 



(dY) = \A*A\^ rn - 1 '>/ 2+ \dX). 
Lemma 2.2. Let S £ ^3^. T/ien ignoring the sign, ifY = S _1 

(dY) = ISI-^-^-^dS). 



(13) 



(14) 



We end this section with a some general results, which are useful in a variety of situations, 
which enable us to tra nsform the density function of a matrix X £ ^ m to the density function 
of its eigenvalues, see iDfaz-Garcfa and Gutierrez- J aimed ( 2009| ). 

Lemma 2.3. Let X £ be a random matrix with density function / X (X). Then the joint 
density function of the eigenvalues Ai, . . . , X rn of X is 



7T 



m 2 p/2+g m 



r^[m/?/2] 



lib n 

Y[(Xi-X S ) p / /(HLH*)(dH) 



(15) 



where L = diag(Ai, . . . , X m ), Ai > • • • > A m > 0, (dH) is the normalised Haar measure and 



Q = 



0, (3=1 

-m, (3 = 2 

-2m, (3 = 4 

-4m, (3 = 8 



3 Riesz distributions 



Altern atively to find the Riesz distribution from the Riesz measure, lHassairi and Lai mi 



(2001), m this section we establish the two versions of the Riesz distribution via the next 



result. For this purpose, we utilise the complexification = 6^ + i&m of 6^. That 
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is, consist of all matrices X G (3 ff ) mxm of the form Z = X + iY, with X, Y G 6^. 

We refer to X = Re(Z) and Y = Im(Z) as the real and imaginary parts of Z, respectively. 
The generalised ri ght half-plane = %L + i&(L in consists of all Z G 6„ e such that 
Re(Z) G qp^, see (|Gross and Richardsl . Il987i p. 801). 

Lemma 3.1. Let X G and k = (fci, &2, . . . , ^m)> &i > &2 > ■ • ■ > > 0. TTiera 
1. for Re(a) > (m - l)/3/2 - fc m , 



ag<p1 



etr{-E- 1 A}|A| -( m - 1 ^/ 2 - 1 ? K (A)(dA) 

= r^[o,K]|sr g)C (s) (16) 

£ for Re(a) > (to - l)/3/2 + fc 1} 

/ etr{-S- 1 A}|A| a -( m - 1 )' 3 / 2 - 1 g K (A- 1 )(dA) 
./A£<p£» 

_ rf>,- K ]|sr 



(17) 



Proof. Let X G and in integrals (Q~6|) and (fT7| make the change of variable A = 
E 1 / 2 WS 1/2 , where S 1 / 2 G *^ denotes the square root of S, such that (S 1 / 2 ) 2 = £. 
By LemmaO (<2A) = \T,\^ m -^/ 2+1 {dW), then integrals becomes 

1. 

/ etr{-W}|W| a - (m - 1) ' 9/2 - 1 q K (S 1/2 WS 1/2 )|S| a (dW). 



2. and 



f etr{-W}|W| a -(" i - 1 ^/ 2 ^ 1 g K (S- 1 / 2 W- 1 S- 1 / 2 )|S| Q (dW). 



The desired results now is follows from ([1]) and ([9]), and (|TT|) and (jTUJ) , respectively. □ 

Hence, as consequence of Lemma 13.11 now we can propose the two definitions of Ricsz 
distribution. 

Definition 3.1. Let £ G and k = (fci, fc 2 , . . . , k m ), k\ > k 2 > ■ ■ ■ > k m > 0. 
1. Then it said that X has a Riesz distribution of type I if its density function is 
1 



■ etr{-S~ 1 X}|X| Q - (m ~ 1) ^ /2 - 1 o K (X)( ( iX) (18) 



for Re(a) > (m - l)/3/2 - k m . 
2. Then it said that X has a Riesz distribution of type II if its density function is 

gw(S) etr{-S- 1 X}|X| a -(™- 1 ) /3 / 2 - 1 g K (X- 1 )(dX) (19) 

r£j[a, -«]|S| a 

for Re(a) > (m - l)/3/2 + fc x . 
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Note that, the matrix multivariate gamma distribution is a particular example of the 
Riesz distribution. Furthermore, if k = (0, 0, . . . , 0) in two densities in Definition 13.11 the 
matrix multivariate gamma distribution is obtained. In addition observe that, in this last 
case if X is defined as Y = 2X the Wishart distribution is gotten. 

In next result is established their characteristic functions of Riesz distributions. 

Theorem 3.1. Let X S & m and n = (k\, k%, . . . , k m ), k\ > k% > • • • > k m > 0. 

1. Then if X has a Riesz distribution of type I its characteristic function is 

\I m -iET|-«?„ (7l m -iS^TE 1 / 2 )^ (20) 



for Re(o) > (to - l)/3/2 - t m . 
2. Then if X has a Riesz distribution of type II its characteristic function is 

IT '" ;STI " ■i m -iiiT\- a q J(i m -ii:^Tj:^)) (2i) 



g K ((l m -iSV2 T SV2) x ) 

for Re(a) > (m - 1)^/2 + fci. 
Proof. 1. The characteristic function is 
1 



/ etr{-(E- 1 -iT)X}|Xr-( m - 1)/ '/ 2 - 1 g«pC)(dX). 



r&[a,/s]|E|«g*(E) y Xe v5 

From (I16|) the characteristic function is 

III- 1 - iTl'^iCS- 1 - iT)' 1 ) = l^-iTSI-y^S-^zT)- 1 ) 
|S|°g K (S) g K (E) 

The final expression is obtained from (JTDJ . The results stated in 2 is obtained in analogously. 

□ 

Finally, it is found the joint distributions of the eigenvalues for random matrices Riesz 
type I and II. 

Theorem 3.2. Let X 6 3?^ and k= (fc 1; k 2 , ■ ■ ■ , k m ), k\ > k 2 > ■ ■ ■ > k m > 0. 

1. Let Ai,...,A m , Ax > • • • > A m > be the eigenvalues o/X. TTien if X /ias a Riesz 
distribution of type I, the joint density of Ai, ... , A m is 



7r m 2 /3/2+ e 

r^[m/3/2]r^[a,K]|S|«g K (I]) 



\ [ xl +k ^ m - 1 ^' 2 - 1 F^{-^- 1 ,'L), (22) 



x 

i=l 

where L = diag(Ai, . . . , A m ) and Re(a) > (to — l)/3/2 — fe m . 

2. Lei #i,...,<5 m , <5i > ••• > 5 m > be the eigenvalues o/X. TTien if X /ias a Riesz 
distribution of type II, the joint density of their eigenvalues is 



m 2 /3/2+e 



— — fi^i-^y 

In „HV.|a 11 V ^ 



r^[TO/3/2]r^[a, K ]|S|' . . 
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(23) 



;.=i 



where D = diag(#i, . . . ,S m ) and Re(a) > (m — l)/3/2 + k\. 
And where 

[ etr{ AHBH* } (dH.) = [ Q F$ ( AHBH* )(dH) 

= „f W ^(A,B), 

see \Diaz- Garcia and Gutierrez- J dimei 1 200i . Eq. p. 18) and Jamek {196 >A . Eq. 60, 
p.483). 



Proof. 1. From Lemma 12.31 

7r m 2 /3/2+ e 



r&[m/3/2]r£,[a,K]|E|»g K (E; 



-1/3 



II( A '-^)" 



Therefore, by @ and ©, 



r^[m/3/2]r^[a, K ]|S|^ K (S)^y l\ 



/ etr{-S _1 HLH*}(dH). 

The desired result is obtained from iDiaz-Garcia and Gutierrez- Jaimed (|2009i Eq. 4.9, p. 
18), see also James! ( 1964 . Eq. 60, p.483), observing that 

f ctrl-E^HLlTHdH) = f o^f-S^HLH^rfH) 



(m) 

0F W^ ( _ S -1 )L) . 



2. is proved similarly. 



□ 



Again, observe that when k = (0, . . . , 0), the corresponding joint densities of the eigenval- 
ues of the matrix multivariate gamma and Wishart distributions are obtained as particular 
cases of (|22|) and (|23|) . indistinctly. 



4 Conclusions 

Finally, note that the real dimension of real normed division algebras can be expressed a s 



powers of 2, j3 = 2" for n = 0, 1, 2, 3. On the other hand, as it can be check in iKabd (|1984l ). 
the results obtained in this work can be extended to the hypercomplex cases; that is, for 
complex, bicomplex, biquaternion and bioctonion (or sedenionic) algebras, which of course 
are not division algebras (except the complex algebra), but are Jordan algebras, and all their 
isomorphic algebras. Note, also, that hypercomplex algebras are obtained by replacing the 
real numbers with complex numbers in the construction of real normed division algebras. 
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Thus, the results for hypercomplex algebras are obtained by simply replacing j3 with 2/3 in 
our results. Alternatively, following iKabd ([198J), we can conclude that, our results are true 
for '2 Tl -ions', n = 0,1,2,3,4,5, emphasising that only for n — 0,1,2,3 the corresponding 
algebras are real normed division algebras. 

Here we prefer our notation used unlike the proposal in Hassairi and Lai mi ( 200ll) . be- 
cause we think that is clearer to a statistici an. Note, however, that any of the results 
obtained can be written easily as proposed bv lHassairi and Laimi ( 2001 ) in the context of 
Euclidean simple Jordan algebras, simply applying ([5]). For example in this case (I2U1) can 
be written equivalently as: 



lK+a (jl 



■ s i/2 TS i/: 
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